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FOREWORD 


Project Whirlwind makes use of the facilities of tho Digital Computer 
Laboratory. The principal objective of the Project is the application of an 
electronic digital computer of large capacity and very high speed (Whirlwind I) 
to problems in mathematics, science, engineering, simulation and control. 

Whirlwind 1 is of the high-speed electronic digital type, in which 
quantities are represented as discrete numbers, and complex problems are solved 
by the repeated use of fundamental arithmetic and logical (i.e., control or 
selection) operations. Computations are executed by fractional-microsecond 
pulses in electronic circuits, of which the principal ones are (1) the flip- 
flop, a circuit containing two vacuum tubes so connected that one tube or the 
other is conducting, but not both; (2) the gate or coincidence circuit; (3) 
the magnetic-core memory, in which binary digits are stored as one of two di¬ 
rections of magnetic flux within ferro-magnetic cores. 

Whirlwind I uses numbers of 16 binary digits (equivalent to about 5 
decimal digits). This length was selected to limit the machine to a practical 
size, but it permits the computation of many simulation problems. Calculations 
requiring greater number length are handled by the use of multiple-length numbers. 
Rapid-access magnetic-core memory has a capacity of 32,768 binary digits. Pre¬ 
sent speed of the computer is 40,000 single-address operations per second, equl- 
valent to about 20,000 multiplications per second. 


1. REVIEW AND PROBLEM INDEX 


This report covers the specific period of March 18, 1957 to May 26, 1957. 
During this time, 60 problems made use of 347.32 hours of the 404.40 hours of Whirl¬ 
wind computer time allocated to the Scientific and Engineering Computations (Sfc.EC) 
Group. Of the 404.40 allocated hours of computer time, 3.3% was down time because 
of computer malfunctions. 43.73 hours of the allocated time were used for terminal 
equipment testing and calibration, demonstrations and various inter-run operations 
not logged to specific problems. 

The 60 problems run during this quarter cover some 18 fields of applica¬ 
tion. The results of 19 of the problems have been or will be included in academic 
theses. In those 19 problems, there are represented 12 Doctoral theses, 5 Master's, 
3 Nuvil Engineer's and 1 Bachelor's. Nineteen of the problems have originated from 
research projects sponsored at MIT by the Office of Naval Research. 

Two tables are provided as an index to the problems for which progress 
reports have been submitted. In the first table, the problems are arranged ac¬ 
cording to the field of application and the source and amount of time used on 
WWI are given. In Table II, the problems are listed according to the principal 
mathematical problem Involved in each. In each table, the letter after the pro¬ 
blem number indicates whether the problem is for academic credit and whether it 
is sponsored. The code is explained in section 2.1, Introduction. 


4 


5 





APPROVED FOR PUBLIC 


PROBLEM INDrfX 


ru id 

Aoroelastlc k 
Structural nest 
Laboratory 

Aar• maut Iral 
Eng t ueor t ng 

Aerophyaics 
Hesearch Group 

Chemistry 

Civil 

Bnyinc*rto( 

Electrical 

Engineering 

Geology 


Hydrodynamics 

Laboratory 

Lincoln 

Laboratory 


Mathematics 


Mechanical 

Engineering 


Motuurulogy 

Naval Supersonic 
Laboratory 

Nuclear 

engineering 

Physics 


Servomechanlama 
Laboratory 



I* rob las 

Minutes of 

Supervisor or 

UascrlptIon 

Nils bar 

WWI Tlss 

Progra—er 

Temperature and atreaa rosponse 

400 C. 

115.*» 

J.C. Lor la 


Horl/.ontal stabilizer analysis 

Blast response of rotor blades 


257 C. 
364 C. 

853.9 

177.2 

E.CriSolone 
E.A. »ltScr 

Extraction of stability derivatives fros flight teat data 

An analytical study of bluff bosb trajectories 

317 C. 
409 C. 

781 .6 
243.0 

L.Maeaola 
F.G Kafir 

Frequency apectrua of sAgnosiua 


4ON C. 

26.5 

L.Slutsky 

Koaponse of a Multistory building undo! 
Dynamic response of shear walls 

dynastic loading 

203 C. 
360 B. 

633.0 

216.9 

H. Gray 

S. Ho 

Switching circuits 


380 B. 

104.3 

C.Roth. Jr. 


MIT anissic project 

Frmrler synthesis for crystal structures 

106 C. 
261 C. 

140.3 

71 .0 

S Sinpsc-n 
M.J.Buerger 

Stokns p.rticle velocities 

J83 C. 

2.3 

T. Marlow 


Eigenvalue prcbles for propagation of elec. rasuiguut it. saves 

193 

L. 

1000. J 

H B. Haight 

Tropospheric propagstiun 

300 

L. 

.8 

*. Mas on 

Error analyals 

312 

L. 

290.4 

L. Peterson 

Prediction analysis 

327 

L. 

700.8 

L.Peterson 

Coverage analysis 

377 

L. 

477.0 

L.Peterson 

Phase error calculation 

395 

L. 

68.6 

D.MacLellon 

Nonlinear second ordor differential equations in the theory 

337 

N. 

81 .2 

H.Wuinitschke 

of elastic shells 

Laminar boundary layer of a steady. compressible floe in 

199 

». 

6.8 

T.Y.Toons 

the ontranee region of a tuoo 

Killing bearings 

293 

C. 

23.0 

A .Shanh.it> 

Growth of futtgue cracks 

3HI 

B.N. 

46.3 

J B.Walsh 

Supersonic flow of air in a tube 

.189 

U. 

134.0 

J.Kadblll 

Statistical and dynamic methods in forecasting 

341 

c. 

547.3 

E. Kelley 

Free convection past a vertical plate 

386 

c. 

153.3 

M. Fins ton 

Diffusion boundary layer 

407 

c. 

33.6 

J . K. Da ron 

Fuel composition m nuclear reactors 

405 

ft. 

15N.9 

M. Cohen 

Theory of neutron reactions 

24 a 

N. 

356.5 

H. Feshbnch 

APB as applied to face- and body-centered iron 

253 

N. 

390.1 

J. Wood 

Analyals of air shower data 

273 

N. 

1235.1 

F. Scherb 

Multiple scattering of waves from a spatial array of 

274 

N. 

71 .3 

P.M Morse 

apnorical scattarers 

Energy levels of dlatoalc hydrldea LlH 

278 

N. 

2384.9 

A, Karo 

Atomli wave functions 

28H 

N 

1222.2 

K. Watson 

Polarizability effects in stuns and aulwculea 

290 

N 

1118 1 

L.C.Allen 

Pure and impure potassium chloride crystal 

309 

B.N. 

300.7 

L P.Howland 

Data reduction 

126 

C. 

794.2 

D.T.Roas 

Temperature distribution in aircraft generators 

J88 

0. 

211.1 

R Moroney 
D.T.Roaa 

Autumn tic programming for numerically controlled machine 
tools 

394 

c. 

517.7 


6 


Table I Current Probleaa Arranged According to rield ot Application 


EASE. CASE 06-1104 


PROBLEM INDEX 

Problem 

Matbssatlcal Problem m-ocedure Wambrnr 

Diagonal laatton HI N. 

Diagonal laaI loo 1M » 

8cha Id 1 process , diagonal lastlonllVO N 

trout 313 L. 

trout 917 C 

Successive approaisatIon. 341 C. 

Croat 

trout, it oration JM C. 

Direct 377 L. 

Met luahey-guine reduction 310 It 

Direr! 3*3 C. 


3. Ordinary differential equation* 

System 

Wave equation 

Different aeta of differential equattooa 
Two simultaneous. non-linear, second order equations 
Seta of simultaneous second order differential squat 1 
Pair of siaultaneoua non-linear equation* 

Hon-1 inter differential equations 
Different aeta of differential equation* 

3 Partial differential equal Iona 

Second order equation Finite difference 400 C. 

4. Integration 

Integration* 

Fresnel integral 

Integration* 

Integration* 

Integration 
Integration 
Integration 
Integration 

3. Transcendental equations 

Curve fitting Lsaat squares 373 N. 

8. Canplem algebra 

Cobles roots and function evaluation Iteration 1D3 L. 

7. Data reduction 

Data reduction Polyncsial fitting, etc. 138 C 

8. Fourier tirloi 

Four lei ayatbesis Direct evaluation 3*1 C. 

9 atattaitcs 

Muittpls lisa aeries Fredlctlos by llnssr operators 10g C. 

10. latogral equations 

Voltarra type Trapetoldal 331 3.N. 


Simpson*a rule 270 n. 

Conversion power aeries 300 L. 

with a cospies argusent 
Barnet and Coulaon expanalon 309 8 

Oauas quadrature 313 L. 

Slapaoo's rule 3M L 

Slap*or'* rule 400 8. 

Kut ta-Gl11 407 C. 

Kutta-GlU 40* C 


Gill 1 D® n 

Ullne predictor-corrector *4ft f 

forsula • 

Step by step 337 C. 

Pd*»f eerie# 337 A.B 

Backeard difference equation# 360 C 

Kutta-0111 306 C. 

Finite difference 300 0 

Kutte-Olll 300 X 


i Salt la algebra end equations 
eigenvalue* 
eigenvalues 

Ortbutuirsal last Ion . eigenvalues 
Inversion 
Linear equations 
Diagonal teat too . inversion 

1 liver*ton. elgnevaluea 
Arvtha*tic operations 
Boolean algebra. Boolean Matrices 
Arltusstlc operations 


Table 11 Current Probleea Arranged According to the Matheaatlce Involved 
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WHlHLVflND CODING AND APPLICATIONS 


2. WHIRLWIND CODING AND APPLICATIONS 


2.1 Introduction 

Progress reports as submitted by the various programmers are presented 
In numerical order In Section 2.2. Letters have been added to the problem numbers 
to Indicate whether the problem Is for academic credit and whether It Is sponsored. 
The letters have the following significance. 

A Implies the problem Is NOT for academic credit, is UNsponsored. 

B implies the problem IS for academic credit, is (^sponsored. 

C implies the problem is NOT for academic credit, IS sponsored. 

D implies the problem IS for academic credit, IS sponsored. 

N implies the problem is sponsored by the Office of Naval Research. 

L implies the problem is sponsored by Lincoln Laboratory. 

The absence of a letter indicates that the problem originated within the S&EC 
Group. 


2.2 Progress Reports 
106 C. MIT SEISMIC PROJECT 

During the past quarter, we have completed the testing of the station- 
arity hypothesis on seismic traces and well logs, and also the construction of 
theoretical random media with controlled scattering properties. Beyond this, we 
plan to do a few miscellaneous computations which may be necessary for the Protect's 
final report. 

S, Simpson 

Geology and Geophysics 

126 D. DATA REDUCTION 

Problem 126 Is a very large data-reductlon program for use in the Servo¬ 
mechanisms Laboratory. The overall problem Is composed of many component sections 
which have been developed separately and are now being combined Into complete pro¬ 
totype programs. Descriptions of the various component sections have appeared in 
past quarterly reports. After the development and testing of the prototype Whirl¬ 
wind programs is completed, the programs will be re-coded for other, commercially 
available, large scale computers, (probably the ERA 1103, IBM 701 and IBM 704 com¬ 
puters), for use by Interested agencies for actual data reduction at other loca¬ 
tions. The programs are currently being developed by Douglas T. Ross, David P. 
McAvinn, Benson H. Schoff and Dorothy A. Thompson, Servomechanisms Laboratory 


staff members with the assistance of John F. Walsh. Tills work Is sponsored by the 
Air Force Weapons Guidance Laboratory through D8A Project 7068. 

Tho nature of the problem requires extreme automat icily and efficiency 
in the actual running of the program, but also requires the presence of human 
operators in the computation loop for the purpose of decision making and program 
modification. For this reason oxtenslve use is made of output oscilloscopes so 
that the computer can communicate with tho human, and manual intervention re¬ 
gisters so that the human can communicate with the computer in terms of broad 
ldeus, while the computer Is running, and have tho computer program translate 
these ideas into tho detailed steps necessary for program modification to con¬ 
form to the human operator's decision. The progiom which does this translation 
und modification is called the Manual Intervention Progran (MIV). The most re¬ 
cent version of the prototype data-reduction program is called the Basic Evaluation 
Program. 

During the past quarter, work has continued on tho development of SLURP 
(Servo Lab Utility Routine Program). As described in the last quarterly report, 
the SLURP program operates like a specially designed simulated computer with un¬ 
usual control features. Work has continued on the latest version of the logging 
program which creates a logical rocord of all manual intervention actions and 
displays during u computer run on magnetic tape. Further additions to the Editor 
Generator program have also bean made and these programs are now in the final 
testing stage. A new Input Translation Program for the SLURP computer, which in¬ 
corporates increased flexibility for making changes in SLURP programs, has also 
been completed. 

Most efforts during this quarter hove gone into the preparation of pro¬ 
grams for data reduction on the Univac Scientific 1103 computer at the Air Force 
Armament Center, Eglin Air Force Base. The Whirlwind 1-1103 Inpu* Translation 
Program developed under Problem 256 for this work has received final modifications 
and extensive translations have taken place during this quarter. It is expected 
that this work will be completed within the next month or so. 

Because of the extreme dependence of the work of Problem 126 on manual 
intervention and oscilloscope output, it is planned that continued use of the 
Whirlwind computer will be made during the next year. Where appropriate, portions 
of the work will be transferred to the IBM 704 computer at the Computation Center. 
But until equivalent specialized input-output equipment is available on that com¬ 
puter, tho basic research in computer applications will continue on the Whirlwind 
computer. 

D. T. Ross 

Servomechanisms Laboratory 

133 L. EIGENVALUE PROBLEM FOR PROPAGATION OF E.M WAVES 

Calculation of eigenvalues, eigenfunctions and mode sums have continued 
for tho bilinear model at 410 and 3000 Me. The 410 Me computations using the pro- 
sent program have been completed; tho 3000 Me computations will be extended to 
higher mode numbers. Similar calculations are being carried out at 50 Me using 
an inverse-square model. 
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Programmer* working on thl* problem are H. B. Dwight and K. M. King. 


WHIRLWIND CODING AND APPLICATIONS 


257 C. HORIZONTAL STABILIZER ANALYSIS 


K. M. Ring 
Lincoln Laboratory 


199 N, 


LAMINAR BOUNDARY LAYER OF A STEADY 
REGION OF A TUBE 


COMPRESSIBLE FLOW IN THE ENTRANCE 




... . -wnere temperature dependenc 

the viscosity and thermal conductivity ol air Is taken Into consideration 
solutions are to be obtained for different entrance .Mach numbers. 


Of 

Mo ro 


203 C. 


T. Y. Toonj 
Mecnanlcal Engineering 


RESPONSE OF A MULTI-STORY FRAME BUILDING UNDER DYNAMIC LOADING 


IXiring the past quarter 
bugging of a new general program, 

cess . 


a five-story program was completed 
using expanded high-speed storage, 


The de¬ 
ls now In pro- 


AddIt tonal production runs have been satisfactorily completed during the 
last quarter. Hopefully, all production runs will be completed by the end of the 
next quarter. 

Programmers working on this problem are N. P. Hobbs and E. Crlsclone. 

E. Crlsclone 

Aeroelastlc und Structures 
Research Laboratory 


261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 

During this quartor we have completed our 3-dimensional transform compu¬ 
tations for the crystal Rhodonite and are presontly contouring the results. Unless 
tho contoured results show behavior which makes us susplcloua of the computations, 
our problem is completed. If there Is trouble, a few computations may have to be 
rerun. 


245 N. THEORY OF NEUTRON REACTIONS 


R. G. Gray 
Civil Engineering 


S. Simpson 

Geology and Geophysics 


273 N. ANALYSIS OF AIR SHOWER DATA 


— • * “b 




* CI 1UU3 


.=.“i 

pleto set of cross-soctIons and angular attributions 'tf , ° n “ com 

been started. ‘ distributions for these values has also 


253 N. 


Elizabeth Campbell 
Laboratory for Nuclear Science 


APW AS APPLIED TO FACE- AND BODY-CENTERED IRON 


and matrix element, for the fac.-emt^'p^r^T'“‘‘ l “ n °‘‘ ri ’‘' 1 coi "‘' letl ° 
Tho computation of tho matrix elements f or P h h f 1,011 " UC nOW bein2 computed 
zation of the secular detain™™^contatnlnu TS”***"* iron “*« diagonal 
face- and body-centered phases still must b^p£t“ * BlenentS for b ° th 1 


Programmers working on this problem are J. H 


■food and 11 , Saffren, 


Most of tho computer time used during the past quarter was devoted to the 
reduction of air shower data. A Monte Carlo analysis of 25 artificial showers is 
in progress for the purpose of determining the resolution of the experiment. As 
yet. we have not started working on the 704 version of the program. 

A paper on the air shower expeiiment will be published soon In Nature . 

The title is An Experiment on Air Showers Produced by High Energy Cosmic Rays . 

Programmers working on this problem are G. W. Clark and F. Scherb. 

F. Scherb 
Physics 

274 N. MULTIPLE SCATTERING OF WAVES FROM A SPATIAL ARRAY OF SPHERICAL SCATTERERS 

2 

Fitting formulae, F(k , K) have been obtained for three additional sums: 
-2ITand -2TT 2 k z Mf e and -2TT 2 k z Mi\j. Each of these covers all seven direc¬ 
tions of K. There now remain five more sums to fit and theae are now being pre¬ 
pared to run on Whirlwind. 


J. H. Wood 

Solid State and Molecular 
Theory Group 


M. Kanakashian 

Laboratory for Nuclear Science 
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278 N. ENERGY LEVELS OF DIATOMIC HYDRIDES LiH 

tl 1 , l! *, ten ’ 1Ve “ belng ■ orte ot “>• “iglt.l computer facilities in ev.l- 

ting all Of the overlap, one-electron and two-electron Integral* which will be 

: “ e0, * u “ the «° rk P- b »« »hich ha. been outlined in pr.vlo^ 

battle h P“ r P°* e including IT orbitals on both centers Is to provide for a 

obta med°before “ DgUlar correl * tion in the molecular wave function than has been 

the WhlrlwIndT*^" ° f the lntegr “ ls is Proceeding for the purpose of utilizing 
the Whirlwind Roothaan programs written by Heckler and Nesbet to determine the 

TtZc tL S lT etTy r Ual “- f ° Und aS Un ° ar conb lnation, of the appropriate 

mlnant representlna th ' ll ? 1Ve th ° l0W08t ° n ° rgy *° r the 8lngle S1,ter doter - 
cMcula^L * r ° Und * tat0 conf ^ration. The SCF HO ground state 

soie ihysrc.r“eaMtv U to th!" 16 ? T “ U " U ° d Cont interaction to give 
tances and r molecular wave function at larger lnternuclear dls- 

, nsequently, to the electronic energy and the dipole moment. 

Programmers working on this problem are A. R. Olson and A. H. garo. 

A. M. Karo 

Solid State and Molecular 
Theory Group 

288 N. ATOMIC WAVE FUNCTIONS 


lowing: 


Atomic calculations with the Roothaan procedure were made for the fol- 


1 ) 

2 ) 

3) 

4) 

5) 


0 with various numbers of basis functions 
(45)* (3d)C configuration of Fe. 

45(3d) n configuration of Fe. 

(3d)“ configuration of Fe. 

(45) (3d)8 configuration of Nl. 


290 N. 


R. Watson 

Solid State and Molecular 
Theory Group 


POLARIZABILITY EFFECTS IN ATOMS AND MOLECULES 


our atomic and molecular problems P ° r For W the r di n t >W 8ngaged ln P r °duction runs on 
we are carrying out a dete™“at!on of the b^^d °" hydr ° een “ UOrlde “ Ol0cul « 
quadrupole moment as a function of internu l 108 ° nergy 1 dlpole moment and 
seeking a Hartree-Foch solut^ o! addUi ° n ' " 

adding basis functions of higher angular momentum^nd f y 8ucco8sl '' el y 

as f orbitals. Besides these calculations we have also ^ *° haVe 8 ° ne 39 “ r 

termlne Hartree-Foch solutions for the free fluorine r at ‘«»Pting to de¬ 

atom. The Roothaan scheme has been used here (as in tie"; f T” lon and neon 
and an indication of the successful convergence of » probiem above > 

recent runs differ from each other by only 03 e ! I that ° Ur " 08t 

proximately 3000 e.v. * J . J e.v. out of a total energy of ap- 


WHIRLWIND CODING AND APPLICATIONS 

Those problems and others developed for electronic structure problems 
have utilised Whirlwind in an efficient and thoroughgoing way. It is also true 
that digital computers and Whirlwind especially have revolutionized the calcula¬ 
tion of electronic wave functions for atoms and simple molecules. In transferring 
to the IBM 704 one can gain the advantage of standardization and it is hoped that 
it will be possible to Interest people in translating and improving programs now 
ln existence. Since this is the final Whirlwind I Progress Report and we have 
not yet finished all of our molecular calculations, wc will continue to report on 
these problems in the Solid State and Molecular Theory Group Quarterly Reports. 

Finally, we should mention a computation underwuy on Problem Number 412. 
No machine time has yet been used but an energy band calculation for Potassium is 
being actively worked on. The calculation will be carried out by the Augmented 
Plane Wave method and will largely use programs built up by M. M. Saffren and J. H 
Wood (these have been discussed previously in Whirlwind Progress Reports). There 
is considerable interest in evaluating the accuracy with which this mathematical 
model represents the physical situation and in particular a variety of assumed one 
electron potentials are to be tested. 

Results of this work will also be given in the Solid State and Molecular 
Theory Group Progress Reports. 


L. C. Allen 

Solid State and Molecular 
Theory Group 


293 C. ROLLING BEARINGS 

During the past quarter, some short runs were made using new input 
data. No future computations are contemplated at present. 

A Shashaty 

Mechanical Engineering 


300 L TROPOSPHERIC PROPAGATION 

An attempt is being made to compute some of the quantities required for 
a study of the bearing of properties of the geomagnetic field upon radio propaga- 
tlon. Programmers working on this problem are P. A. Duffy and J. N. Skenlan. 

P. A. Duffy 
Lincoln Laboratory 

309 BjN. PURE AND IMPURE POTASSIUM CHLORIDE CRYSTAL 

Details of the charge density of KC1 in an x,y,0 plane containing nuclei 
nave been obtained. Now the resulting contours of the charge density have been 
drawn and comparisons with recent experimental data have been made. Information 
on a few other minor points may still be sought on Jfhlrlwlnd as long as the ma¬ 
chine is available. 


L. P. Howland 
Physics 
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312 L. KKHOR ANALYSIS 

This problem, with its various modifications, has boen used for produc- 
tlon runs during the past quarter. There .111 probably be more production runs 

son and L! U peter.on nC ° ln ,lm °' PrOgr “" 0rB *° rkln « thls E. Hutch.- 

L. Peterson 
Lincoln Laboratory 

317 C. EXTRACTION OK STABILITY DERIVATIVES EHOH PLIGHT TEST DATA 

lateral stahiM? P“ at 1 uarter ‘»ls project has continued the investigation of the 

v‘tlve mlih^ r “ 8 eXlr “ ction P~ bl «. “Win* Shinbrofs method, the deri¬ 

vative method and a new statistical method developed by this project Investlea- 

n *f° hrr ,,C , t °* ? U8tS accuracy for longlLinal 

srs?-■* *■ -•*“ >-• «- 

ino investIgatlun, non tn progress. Kesults ntll be potjltnhso" 1 Uw'ho rl,“£ letin * 
Wright Air Development Center Technical Notes. * 

u so Sn working on this problem are M. E. Hoult, L. L Mazzola 

M. N. Springer, J. F. White and L. E, Wilkie. ' 


M. N, Springer 
Aerophvaics Research Group 


PREDICTION ANALYSIS 


3) How various Input parameters affect the functioning of the system. 

I 1 " the*future? Vei,t ** at * ‘ f8 * "° r8 '—«.rl.tlc. . there win be some testing done 

= ,n. n.,H the . L n„0 S.""” “ “r“" 

The program concerned with the effects of the 

“ unB ■*»* 

und.r.t„d: n r^ c “::n:rth7d:L s ;^sS^^ d of ^ 
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WHIRLWIND CODING AND APPLICATIONS 

mentioned in Progress Report No. *10). This now program Is also interesting In that 
it has been written to handle a much more general sot of conditions than has been 
considered in previous programs. 

Programmers working on this problem are S. Hutcheson, L. Peterson nnd 
P. Willmuiin. 

L. Peterson 

Lincoln Laboratory 

337 N. NONLINEAR SECOND ORDER DIFFERENTIAL EQUATIONS IN THE THEORY OF ELASTIC 
SHELLS 

Computations for problem 337 have been completed this term. The results 
have been written up for a Technical Report, submitted to the Office of Statistical 
Services. 

II. Wuinitschke 

Mathematics Department 

341 C. STATISTICAL AND DYNAMIC METHODS IN FORECASTING 

During this past quarter, work has been completed on a project dealing 
with a statistical-dynamic approach to numerical weather prediction. The results 
wore not especially encouraging, however, the changes necessary in order to raako 
the method successful are believed to be known. 

A short investigation of the possibility of forecasting mean monthly 
temperatures at « single station (Washington. D.C.) for a given year, from the 
empirical orthogonal functions of the mean monthly temperatures at the same sta¬ 
tion for the preceding year, was also conducted. Preliminary results are quite 
promising and may be of considerable Interest to those engaged In long-range 
forecasting. 

Programs for testing developed linear forecast methods on a set of inde¬ 
pendent data were completed and run during the lust quarter. This brought to an 
end this particular phase of studies of sea-level pressure by empirical orthogonal 
fund ions. 

Computations have also been completed in an Investigation of a theoreti¬ 
cal model of the atmosphere using an expansion of the meteorological variables in 
a small number of spherical harmonics. In a final run the system of equations was 
numerically Integrated by a stepwise method over the equivalent of two hundred 
days. 

The bulk of the computations Involved in the studies mentioned above 
have been perfumed on Whirlwind I. The programmers for this problem are: E. 
Kelley, B. Short , W. Sellers, K. Dryon and J. MacDonald. 

E. Kelley 

Meteorology Department 
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WIIRUVIND CODING AND APPLICATIONS 
360 B. DYNAMIC RESPONSE OF SHUAR WALLS 

square) <,ul,rtor ' rosult8 *ero obtained on 3 x 2 grids (22 ' 

Q /. fforts have been made to obtain results on 6 x 4 grids. 


Shui Ho 

Civil Engineering 


3G1 B,N. GROWTH OK FATIGUE CRACKS 


=s •?£=•=■£HS™™A=r=“ 

ys=S535sS«i 

J. B. Yalsli 

Mechanical Engineoring 

334 C. BLAST RESPONSE OF ROTOR BLADES 

rotor has boon'dcuusee^InU Is'readyT *' 0 Shai,OS and fr equencles of the 
being written to c™ ' e S L n ' f °” rUnS ' -re currently 

tbe expressions for the ieneraUzed alrf orL’ p!'^!' Certa , ln C ° efflclo "‘ 8 

x^sinr*r tion uf " t - 

study of blast response of helicopter roun s'wtu'then ^0^04; A 1>aran,etrlC 
Programmers working on this problem are B. A. Vltmer and A. Holland. 

E. A. tfitmer 

Aeroelustic and Structures 
Research Laboratory 

377 L. COVERAGE ANALYSIS 

used in i,eP ° rt NO ' 4,5 ' haB bee " 

useful for the annlysls desired i ,,4,, . 11,0 ro9ultll *t: Rata has been found 

for some time. ■ U WiU continue to be n production run problem 

Programmers working on this problem are E. Hutcheson and L. Peterson. 

L. Peterson 
Lincoln Laboratory 

380 B. SWITCHING CIRCUITS 

" — VrVzti sxz 
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written up in a S M. thesis entitled. Digital-Computer Solution of Switching- 
Circuit Problems . Four Whirlwind programs have been written for solution of 
switching problems: 

1) Analysis of Single-Output Contact Networks. The transmission of 
any contact network with 3 to 42 nodes and 4 to 8 variables can 
bo determined by this program. The network is specified by u 
Boolean matrix and the transmission is in standard sum form. 

2) Analysis of Multiple-Output Contact Networks. Similar to (1) except 
that the transmissions between all pairs of nodes in the network 
are obtained simultaneously. 

3) Synthesis of Non-Serios-Parallel Contact Networks. This program is 
capable of synthesizing a contact network to realize any four- 
variable switching function. The desired transmission is specified 
in standard sum form, and the output is a Boolean matrix. 

4) Determination of Prime Impllcants and Minimum Sums. Given the 
standard sum and optional torms (maximum 512 terms) for any switching 
function of 14 or fewer variables, this program will determine the 
prime impllcants and one or more minimum sums. 

Using the latter program, approximately 100 problems were solved for J. P. Gouyot 
in connection with the development of a now character decoder for a photo-composing 
machine. 

C. II. Roth, Jr. 

Electrical Engineering 

333 C. STOKES PARTICLE VELOCITIES 

Problem 383 was the solution of Stokes equations for the water particle 
velocities in an oscillatory wave. The equations were solved for 24 values of 
wave phase angle and three depths for ‘11 waves. The wave parameters and positions 
were identically those used in an experimental determination of particle velocities. 

The experimental and computed values were plotted to determine the range 
of applicability and adequacy of Stokes relations. Tills report served us u Mister 
of Science thesis and D.S.R. report. 

The solution of Stokes equations by Whirlwind allowed this study to be 
mado. since manual computation would have been impossible from a time standpoint. 
The entire study used less than an hour of machine time with the main routine 
coded in the binary system. 

T. Marlow 

Hydrodynamics Laboratory 

386 C. FREE CONVECTION 

Computations were completed during the present quarter. Exact solutions 
were obtained for 12 representative surface temperature distributions of the form 
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1 + Bx q (for B = 0.1, 1.0 and 10.0 and q ■ -.01 - I - 2 an d - 25) Th» . 

ndloate lh e r.l.tlve effoci, (on skin friction'and heat tranafeof alhanaon 
q to bo insensitive (on order of 0.1%) to changes In I). " fc ° ‘ n 

J. It. Huron 

Naval Supersonic Laboratory 

3«7 C. DETERMINATION OF VELOCITY POTENTIAL 

>egun on T T" C ° raP ‘° U,d f ° r 1,10 dolt “ «*"« work has 

uf the calculations ..a^o ^on^r:^" h"” ““ ^ 

the swept w^“ r on P !he%o; e if U Uce n ss S a h r; he 8, "* 1 ‘ !ht ^ "" “ — add - 

Programmers working on this problem are a. Stnpleford and ». Wonthcrhill. 

H. Stapleford 

Aeroelastic and Structures 
Research Laboratory 

388 D. TEMPERATURE D1STKIDUTION IN AIRCRAFT GENERATORS 

hoped that C cer^In reflioLnrs^'boinc^^urarje^thr 8 ^^ 15 ' “ 

case, this problem will terminate on August 1 . situation, but in any 

K. Moroney 

Servomechanisms Laboratory 

389 D. SUPERSONIC FLOW OF AIR IN A TUBE 

.... p.»>-.....,««, „ „ 

vlty, recalculation of the first aDoroxlMti P ti ° n s. Because of this sensiti- 
for the second approximation have been undertaken •1™^^*°" °* calculatlo " s 
atlon. The tapes for these calculations have 1 on ? oxact proport l' va ' - l- 

wlll require only minor modlficon..ns completed and checked out and 


Mechanical Engineering 

394 c. AUTOMATIC PROGRAMMING FOR NUMERICALLY CONTROLLED MACHINE TOOLS 

lntr«luced T :; e „eTe P r:“rthe l"« ““chine in 1951 

rushed parts. The primary problem associated with (h^tvo re ?“ i, ' e " ,ents lnto «- 
preparation of detailed instructions for m „! ° f mochine .ool, the 

Whirlwind Compute.. Under problem 132 (WhirlwlnI E ' *“* bB °" studled cslng the 
ol subroutines was developed to assist K'l "‘T* “•">«» 33-43) a llbrar 
•roblorn 250 (Whirlwind Progress Reports 41 42) E! I ?" ^ th< “ MIT "“»• d " d « 

the feasibility of applying modern 
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• ii.nnatic programming techniques to machine tool programming was successfully de 
nwnutr.U'd for plane purls bounded by straight line and circular segments Problem 
394 . s concerned with the extension and elaboration of these automatic programming 
techniques to arbitrary surfaces which can be muchined by 3 and ' axis numerical lv 
>n t To 11 od mil ling machine* Tito programs are currently being developed by Dougins 
r. Iuhs , Donald F Clements, Harry E Pople, Jerome it A?nUor. Samuel M. Matsa. 
t?rv one chan isms Laboratory staff members. Tills work is sponsored by the Ur Ma¬ 
teriel Command through U*W Project r»873 

The problem may be viewed as the design of an APT System (Automatically 
Pro.;rammed Tool system) combining the human operator, the machine tool itself and 
m ‘ PT >>mpiiter, which is u combination of a general purpose computer and existing 
machine tool directors The objective is a system of Whirlwind computer programs 
which will simulate an \PT computer with a language which is convenient for human 
expression of machined part designs. A design expressed in this language will be 
automatically translated by the apt computer into machine tool instructions 

Lkiring the past Quarter, a new systematized solution to the problem of 
calculating cutter motions for the *PT 11 system, which allows programming of a 
part by means of arbitrary curves on the surface of the part, has been made. Out¬ 
put routines developed under Problem 132 have also been incorporated to provide 
punched paper control tape for the MIT Numerically Controlled Milling Mhchine at 
the Servomechanisms Laboratory, so that test parts can actually be machined. The 
new method of calculation depends only upon two programs for each type of surface. 
The normal vector program produces a vector perpendicular to the surface of the 
part at a specified point. The directed distance program calculates the distance 
to the surface from an arbitrary point in space looking in an arbitrary direction. 
Normal and directed distance programs for planes and spheres have been in use for 
•ome timi These programs are being combined Into a generalized program which 
will calculate normal and directed distances to planes, spheres, cones, cylinders 
and general quadric surfaces In three dimensions. This program is now under test. 

The programming for the APT III system, which allows programming a part 
by means of entire regions of the surface of the part, is also under test. The 
analysis used in this program also represents a systematized solution to the pro¬ 
blem and depends only upon normal vector and directed distance programs. a number 
of successful tests have been completed and the program Is undergoing further re¬ 
finement. It is planned that, with very slight modification, the programming for 
the APT III system. The oscilloscope output routine which draws noar perspective 
views of the part being made has been used with all programs. 

The work on Problem 394 will transfer to the IBM 704 Computer at the 
Computation Center in the very near future During the week of May 20 a meeting 
was held at M I.T to coordinate the efforts of a number of aircraft plants in a 
Joint programming effort to develop an APT II system for two dimensional parts, 
ulong the lines set forth in these reports. The Servomechanlsms Laboratory, M.I.T.. 
Is to act us coordinator in this Joint effort, and It Is hoped that a preliminary 
system for the IBM 704 can he completed by October 1957. Concurrent efforts, 
closely coordinated, will develop a similar system for the IBM 650. 

D. T. Ross 

Servomechanisms Laboratory 
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395 L PHASE ERROR CALCULATIONS 

Bent “ l “' been fOUnd that the Pr ° ble "’ **■ pr.- 

n has ZTZIT 1 ot WWI tlBe - T —' 


400 C. TEMPERATURE AND STRESS RESPONSE 


D. MacLcllan 
Lincoln Laboratory 


performed ““ ddltl ° nal hlad ° “• Programmed and the run 

ccessruiiy. Currently, three other blades are being programmed. 


405 B. FUEL COMPOSITION IN NUCLEAR REACTORS 


J. C. Loria 

Aeroelastic and Structures 
Research Laboratory 


following e T q h ;.t l ?r‘ Ctl0n rate betW ° en neUtr ° n8 * r ‘ d nucl « 1 «« ^ described by the 


dN 

dt 


rh 


§ N 


( 1 ) 


the cross It* *“ a ,Unctl °" ° f the "«*tron energy. .. „ 

and flux in eq Q) " UC h t ^ t neccssar >’ to d «^n* the effective cross section 


$ '<p (B)d£ 

f-X <r(E)(J) (E)dE 


( 2 ) 


(3) 


the cross lUTon Ir.ss'u^dTvary'“^“rT' 1 ^-“oU.mann distribution, and 
assumptions under certain conditions) then e^ ,1“!™ veloclt y (both reasonabl 

n«n eq. ( 2 ) and (3) can be evaluated exactly 

distribution ^rte C ie e d ar by r °;he°. r f?:ct t s he o r f aSS T Pti0n “ va “ d - Th. flux 

of neutrons from fission energy The crons'^ect Ln^T' 10 ”’ a " d the 8lowln 8 down 
in reactors do not follow the 1/v law. ' Sectlons oi » number of nuclides found 
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if the neutrons were in true thermal equilibrium with the surroundings. Neutron slow¬ 
ing down theory will be used in the range from 6 ev. down to predict the flux. 

From (the upper limit; 6 ©v.) down to £2(5kT, where kT is the peak of the 
thermal Maxwell-Boltzmann distribution) the epithermal flux calculated from slowing down 
theory will be simply added to the thermal flux distribution to give the resultant flux. 

The epithermal flux calculated from slowing down theory approaches infinity 
as the energy approaches zero and at energies of the order of kT is obviously inap¬ 
plicable. Therefore below E2(5kT) the epithermal flux calculated from slowing down 

theory is multiplied by a function Q (see eq. 6e) which reduces the epithermal ilux to 

zero below kT, and provides a smooth transition to the Maxwell-Boltzmann distribution 
in the range from 5kT to kT. Below E 2 the total flux is assumed to be the sum of the 
Maxwell-Boltzmann flux plus this corrected epithermal flux. 


With the flux so defined, as the sum of a Maxwell-Bol tzmann and a slowing 
down contribution, eq. 3 may be written as: 


r E > fi 

C7- (E) 0 mb (E) dE + I CT(E) 0(E) dE 

° $MB Kr $MB 


(4) 


The first term is the contribution from the Maxwell-Boltzmann flux, the second the 
contribution Iron the epithermal flux. From the first term, as compared to the values 
calculated, assuming 1/v cross sections and a Maxwcl1-Boltzmann flux, the effect of 
neglecting the non 1/v cross section characteristic may be determined. From th© se¬ 
cond term may be determined the effect of neglecting the epithermal flux. Both of 
these compoenets of the effective cross sections are numerically evaluated using known 
values of neutron cross sections as a function of energy and the Maxwel1-Boltzmann and 
epithermal flux distributions which are calculated lor the particular reactor condi¬ 
tions under consideration. The flux distributions must be calculated using the fol¬ 
lowing equations. It 1 b a function of the fuel composition and thus will change dur¬ 
ing reactor operation. The primary object of the calculations is to determine the re¬ 
lative magnitude of the contribution to the effective cross section due to the epither¬ 
mal flux, and its variation with fuel complsition. 


Calculation of Epithermal Flux 


0 (E) 

* 


P(E) 


+ r^(E) 


(E) Ci 


where the various terms are defined 


by: 


(5) 


numerical methods. P™Po cd calculations is to solve eq. 


(2) and (3) by 


fi^ov and the upper l^« at about 2 Mev. Between 

r::rr^ ri 

range Is assume to follow „ fax veil Bolt-mann dlstrlb 1°"' ^ flUX in the therm: " 

peraturo Is r ol .ed to account for the ch^r~ IhTd effective tern 

"■■uajuuon to be expected 


F(E) S 


p(E) 


-Op Vjp C 
( -<* 


Cj S 


OC 




fi(E) S p(E) (6a,b) 

f 

(6c,d) 


VrT 
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The computations involved are as follows. The resonance oscnpe probability, 
p(£). must be evaluated as a function of energy. Then the parameter alpha must be 
evaluated. The constant Cj can then be determined. The epithermal flux as a func¬ 
tion of energy can then be obtained from ea. (3). This is then used to evaluate 
the second term in eq. (4), which then gives the required result, for one nuclide. 

It is estimated that 140 to 150 steps will be needed in the numerical in¬ 
tegrations. About nine nuclides will probably be considered, und since both ubsorb- 
tion and fission cross sections are needed in some cases, on the order of 13 x 140 
or 1820 cross section values will be used. 


The contribution to the total effective cross section from the epithermal 
flux may run around fifty percent in some cus-^s, due to large resonances in the 
epithermal energy range. 


IXiring the past quarter, the effective absorption cross sections of 
U-235, U-2J6, U-238. Np-237, Pu-239, Pu-240, Pu-241, Xe-135, and Carbon, as well 
as the effective fission cross sections for U-235, Pu-239, and Pu-241 were calcu¬ 
lated in eight different cases. A sodium-graphite reactor wu; used as the basis 
for calculation, and effects of varying fuel composition, temperature, and macro¬ 
scopic scattering cross section of the system vore considered. 

It was found that neglecting the epithermal flux led to estimates of the 
cross section that were Lower than the true effective cross section, in some cases 
the effective cross section being eight times as large (for Pu-240). 

It was found that the effective cross section was significantly affected 
by variations of ail three parameters mentioned above. Therefore in any calculation 
of fuel composition over extended irradiation times it will uo necessary to recal¬ 
culate effective cross sections at suitable time intervals, on the basis of the cur¬ 
rent fuel composition. This type of calculation, using the methods developed for 
the present problem for calculating the effective cross sections, is being considered 
for future work, using the IBM-704. 

M. Cohen 

Nuclear Engineering 

407 C. DIFFUSION BOUNDARY LAYER 

Mass transfer into a high speed boundary layer results in the simulta 
neous action of momentum, energy, and mass transports and produces considerable 
changes in the thermal behavior of the layer. Previous computations (Problem 207) 
indicated the effects for free flight conditions. Presently the effects to bo 
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encountered under wind tunnel conditions ure t t>o luu ted in conjunction with 
an experimental program into the phenomena. 

The descriptive system of simultaneous equations is of the form 


(momentum) 

(mass) 

(energy) 


(\f")' + ff" - 0 
Cl" + acj ’ a 0 
T" + bj T* 4- b 2 


in which a, , and bo are functions of the Indicated dependent variables The end 
point boundary conditions Involve coupling between cj and X and otherwise constant 
specifications for either the dependent variable or itB derivative. 

The Glll-Kutta routine will be employed for integrations and repeated In¬ 
tegration continued until proper end point conditions art* attained. Successive os 
timates of boundary conditions will be made on the basts of an influence coefficient 
routine supplied to the computer. 

Integration will be over the range U* t) £k.O in increments of A ^ 0.1 

Output will be taken on magnetic tape and oonslus of f, f», f“, ^, ej T, T‘ all 
as functions of tj . 

During the past quarter preliminary computations were made to chuck out 
the program and to evaluate the need for special Integration intervals in regions 
of rapid variations of the parameters. IXiring the next quartor the computations 
will be devoted to obtaining solutions corresponding to wind tunnel test conditions. 

J. R. Haron 

Naval Supersonic Laborutory 


408 C. FREQUENCY SPECTRUM OF MAGNESIUM 

The development of digital computers capable of making the requisite 
numerical calculations has led to u revival of interest in the Born-von Karroan 
theory of lattice dynamics. The vibrational frequency distribution function and 
the bulk thermodynamic properties of several metals crystallizing in the face cen¬ 
tered cubic system have been computed using interatomic force constants derived 
from the adiabatic elastic constants of a metal single crystal or from the tempera 
ture diffuse scattering of X-rays. Since the most sensitive comparison with expert 
moit is provided by the lattice heat capacity at low temperatures, calculations 
Ui,ing the force constants at liquid holium temperatures are of the greatest interest 

Recent low temperature measurements of the elastic constants of magnesium 
have been interpreted in terms of a central force model considering interactions 
between first, second, and third nearest neighbors. It is proposed to solve the 
secular determinant derived from this model at 793 points In the reduced Brillouln 
zone. Since magnesium is a hexagonal close packed metal with two atoms per unit 
cell! the secular determinant is of order f> and 4758 distinct roots will be ob 
tained. It is expected that this root sampling, when combined with the group theo¬ 
retical and topological considerations regarding the existence of singular points 
in the frequency distribution spectrum discussed by Phillips, will provide a good 
picture of the density of vibrational energy levels in magnesium. 
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The computet Iona' problem Involves the generation of the elements of the 
secular determinant at each point in reciprocal space whore solutions are desired 
and the determination of the characteristic values and vectors making use of the 
program developed by F. J. Corbato. 

[Airing the past quarter a satisfactory program was arrived at. The compu¬ 
tations are in progress. A trial spectrum calculated from the roots of the secular 
determinant at 143 points in the first Brlllouin zone have heat capacities in good 
agreement with the experimental data. Programmers working on this problem are E. 
Mack, L. Slutsky and C. W. Garland. 


L. Slutsky 

Chemistry Department 

409 C. AN ANALYTICAL STUDY OF BLUFF BOMB TRAJECTORIES 

A machine solution is to be obtained for a set of 3 differential equa¬ 
tions describing the rigid bodv longitudinal dynamic motion of a bluff-bomb con¬ 
figuration from the instant of release at an initially prescribed altitude, atti¬ 
tude, and speed until it contacts the ground. Non-linearities of aerodynamic 
coefficients with speed, altitude, and angle of attack will be represented ana¬ 
lytically from point to point of the trajectory and the solution will be integrated 
as a function of time with respect to a Euler reference axis system. Approximately 
300 Kutta-Glll solutions of the set below will be required: 

u = -w 0 + C (u,w) gS Q - g sin © 

W 

w a u 0 + C (u,w) fcS Q + g COS (h) 

* w 

9 = C n8 <u ' w) S d2 3 9 + c _ rf (U.w) S d 2 Q w 

21 u 21 u 2 

y y 

+ C (u.w) Sd q 

m — * 


The Greek notation consisting of: 

cC Alpha, lower case angle of attack 

Theta, lower case pitch angle 
© Theta, Capital path angle 

The Roman notation being: 

u Velocity along X-axis 

W Volocity along Z-axis 

c x X-Force aerodynamic coefficient 

C z Z-Force aerodynamic coefficient 

c m Moment coefficient in pitch 






C * Aerodynamic Stability Derivative wtih respect to pitching 

velocity. 

C * Aerodynamic Stability Derivative wtih r«9pect to rate of 

change of angle of attack, 
g Acceleration of gravity. 

S Reference area, 

d Reference diameter. 

W Weight. 

I Pitching Moment of Inertia. 

y 

q Dynamic pressure. 

During the past quarter, the present phase of the problem was completed. 
Future pltns will depend upon the analysis of the present results. 

Programmers working on this problem uro P. G. Kase and T. Carney. 

P. G Kase 

Aerophysica Research Group 
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3.1 Introduction 

analysis, Tt“ll *" C °" PU “‘“°"' n—ricl 

Thu present list of subjects^ dlrecUy relatcTt T *™ 1 CUrr ‘-' ntl > <>«-ed at M.I.T. 
the Academic year of 1957 , r , H { 1 to machln « computation offered during 

/car or 1957 Includes the following: 

Subject Description 


6.25 
6.535 
6.538 
6.5-1 

6.567 

6.568 
2.215 
15.542 
M39 
M411 
M412 


■lachlne Aided Analysis 

Digital Computer Coding 

Electronic Computational Laboratory 

Pulse Data Systems 

Switching Circuits 

Switching Circuits 

Methods of Engineering Analysis 

Management Information Systems 

Methods of Applied Mathematics 

Numerical Analysis 

Numerical Analysis 


Year Instruct or 

•• Arden 
0 Arden 
G Verauh 
0 Wood row 
G Caldwell 
G Caldwell 
G Crandall 
G Gregory 
G Hildebrand 
G Hildebrand 
G Hildebrand 


It is apparent from the above list that 

subjects. However, a number of undergraduate students *” P ™ doffll Graduate A 

on an elective basis. B 1 students may and do take these subjects 


32 Electronic Computational Laborator 


The subject matter offered in k sia .1 

vided the student with an opportunity to ob^af T* ° Spring terB oi 19 ” pro- 
tlon of the Type 650 Magnetic Drum Calculator F ^ kn ° wled * e of the opera- 

in the laboratory operating and using .he BSo'calculator^" 1 SP<!nt “* 1<!a8t 14 hours 

With the basic esHanguage!" Afteethe bl^cTang 6 ^ b> ' flrSt becomln e *««Ul«r 
taught to use the M.I.T. Selective Svo™ ^ FanKua « e was understood, they were 
Pretive Numeric Coding System. On occas ion"cert i FiX ° d a " d Floatln ^ Pol "t Inter- 
one of the less common Interpretive systems ! 1 StUdents found It helpful to use 
Decimnl Interpretive Mneumonic Codlng^ystem and ^ ' For exam P le the MITILAC 

Matrix System were used on a very lilted “.u" ' F1 ° at *'>e Interpretive 

3.3 6_^3 0 ^ nj e j ^ ii o ^^ the SDr ,„ ff T _ 

For a number of years it 

as one of the major undertakings In this .^Wecr“°n P f netice to assl * n a t.r» P a Per 
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the linear scale. The following term paper assignments were performed as Indicated below 

The Feasibility of Monte Carlo Matrix Inversion Using 
Digital Computers. 

Modified Aitken Method of Matrix Inversion. 

Power Limit Transmission with Static Load. 

Sorting on the Dutamatic Computer. 

A Study of Brokerage Accounting. 

M1TSS Subroutines. 

An Attempt at Finding the Characteristic Equation of 
a Complex Matrix. 

Bata Processing in the Brokerage Field. 

Floating-Point Arcsin Subroutine Trace Subroutine. 

Complex Arithmetic Subroutines for Use with M1TSS Floating 
Decimal Point System on the IBM 650 Computer. 

An Integrated Program for the IBM 650 Magnetic Drum Data- 
Processing Machine to Calculate Certain Statistical 
C^iantlties from a Group of Data. 

Subroutines for 18 Digit, Floating Decimal System. 

Motion of a Body in the Earth-Moon System. 

A Digital Computer Solution of Langhar's Velocity Distribution 
Equation. 

A Subroutine for the Calculation of Statistical Variations. 

Feasibility Studies of Automatic Data Processing Equipment 
in Business 

A Study in Market Research. 

An Exploration of a Logical Model to Obtain an Initial Distribution 
in a Stepping Stone Linear Program. 

A Study of Payroll Calculation on a Computer. 

Factorization of a Polynomial (Hitchcock’s Method) 

Double Register, 18-Dlglt, Floating Point Arithmetic Routines for 
the IBM Type 650 Magnetic Drum Data Processing Machine. 

Preliminary Report on Function Display Program for the 704. 

Digital Differential Analyzer Solution of Several Differential 
Equations. 

Market Analysis of Electronic Data Processing Machines. 

Solution of Bessel's Equation by Prodictor Corrector Type 
Formulae on the IBM Type 650 Computer. 

The Solution of Simultaneous Linear Equations on the IBM 704 
Electronic Data Processing Machine by Use of Aitken's 
Elimination Method. 

Routine for the Solution of Any Second Order Differential 
Equation. 

Calculation of the Eigenvalues of a Matrix in Jacobi Form on 
the IBM 650 Calculator. 

Error of Runge-Kutta Fourth-Order Formula Applied to a Bessel 
Equation. 

Programming of Statistical Quantities for Computation on the 650 
Magnetic Drum Calculator. 

650 Program for Solution of Bessel Equation. 

VADFAN. A Program to Perform Variable Degree Factorial Analysis. 

A System of Matrix Arithmetic Subroutines for the IBM 650 Duta- 
Processing Machine. 
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Economical Input to the Machine , i 

^..pantoents Operation on a Pfc.^ , [nfor » al ‘° n re¬ 
solution of the Heat Flow Equation 

COlCU ims n and S^P In8Ur “ nCe PreBlU " S U8ln * th ° 650 with 
Uata -*• —«- —. 
pi'srs: i e xzz™ mgitai — 

Simulation of Consumer Goods Distribution. 

Sorting".^! «or^n r g n on a .M^a U c“ 0 and Wl n “ Err ° r Crlttrl °" 

An Interpretive System for Double Precision Arithmetic. 
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Khode Island, University of 
Tufts University 
Vermont, Univorstiy of 
Wellesley College 
Wesleyan College 
Williams College 


wiiijiuianon umei 


and Associate*! *t ,h«« 7 177 7 of 24 Kesearrh AssUtwi 

cse institutions during the academic year 195(5 - 1957 The«*j 

anS° rr VT Ae aCt,VU llalSOn bt ’ ween lh ° «*« - the Lnl.r and “I’.t^e”. 
at thelr individual Institutions. Appointments of thts type will be 
c year to insure a widespread und dynamic participating program. 


3 1 


Computation Center 


Compton Laboratory and ' ha ’‘accLp^Iei battaliaii C °° Pletlon ot ‘he new Kari T. 

ter in the new building. nation of the large-scale Type 70-1 Corapu- 

modern co»puting C “uIpmenrand n a e brMd“ew^hy»Icar h M h “ US “ ny assetS: staff 

Center represents a cooperative activity invo^vi^u.T 8 ^ 1 * 1C “ Uy 1 thL> Computation 
P-aent. 25 New England Colleges ^d ^versuL"* ' IM C -porat ion'and, “ 


Aocfc iULuieu in me 






-s-www .v«. » v-MH.jjvun LUUOrjlorv . 

second",^ ° f lht ' bas< ' m ‘"’ t - the entire first floor, and part of the 
house t^ ?BM°L r Pl0 " , '“ b0rat0ry - in addition, a two-story anneiTiJ used to 

house the IBM Type 704 Electronic Data Processing Machine (EDPM) and the associated 
Electric Accounting Machine (EAM) equipment. 

The first floor contains adequate space for the headquarters staff, the 
operations staff (analysts.programmers machine operators, etc.), IBM Institutional 
Representatives, New England Unvorslty Research Assistants und Associates, MIT 
Research Assistants and Associates, classroom and seminar room, as well as the 704 
computer. The basement provides space for the EAM machines, the systems research 
aboratory, dark room facilities, the electrical power plant, and the air conditioning 
equipment. The second floor provides space tor the programming research staff, the 
visiting professors, nnd the library and document room. 

All this area has been furnished in a first-class manner to facilitate 
the progress of research at the Center. 




Amherst Collogo 
Bates College 
Bennington Collego 
Boston College 
Boston University 
Bowdoin College 
Brandeia University 
Brown University 
Connecticut, University of 
Dartmouth Colloge 
Harvard Unlvorsity 
Maine, University of 
Massachusetts, University of 
Middlebury College 
Mount Holyoke College 
New Hampshire, University of 
Northeastern University 


The computational.facilities in the Center are supported in large measure 
by the IBM Corporation. Specifically. IBM is providing the 704 computer, the asso¬ 
ciated EAM equipment, and the associated maintenance personnel on a gratis basis. 

The following machine complement is available ln the Center: 


Quantity 


Description 

1 

70*1 

Analytical Control Unit 

1 

711 

Punched Card Reader 

1 

716 

Alphabetic Printor 

1 

721 

Punched Card Recorder 


733 

Magnetic Drum Unit (8192 words) 

1 

736 

Power Frame *1 

2 

737 

Magnetic Core Storage (8192 words) 

1 

740 

CRT Output Recorder 

1 

741 

Power Frame 02 

1 

746 

Power Distribution Unit 

1 

753 

Magnetic Tape Control Unit 

10 

727 

Magnetic Tape Units 

1 

780 

CRT Display Unit 
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Off-Line Equipment 


Quantity Type Description 


1 

714 

Card Reader 

1 

717 

Alphabetic Printer 

1 

722 

Card Punch 

2 

727 

Magnetic Tape Units 

1 

757 

Printer Control Unit 

1 

758 

Punch Control Unit 

1 

759 

Auxiliary 

Card Reader Control Unit 

Machines 

1 

024 

Key Punch 

5 

026 

Key Punches 

3 

056 

Verifiers 

1 

066 

Printing Card Unit )Daia Tran. 

1 

068 

Telephone Signal Unit) & Rec. 

1 

077 

Collator 

1 

082 

Sorter 

1 

407 

Accounting Machine 

1 

519 

Reproducer 

1 

552 

Interpreter 

tXiring the 

Spring Term of 1957 

the 704 Computer and the physical plant 

were not sufficiently 

well along to permit actual use of the machine. As a result 

the Ptudents of 6.538 

weie given guided 

tours through the Center and had an opportunity 

to view the physical 
use the computer for 

machine components, 
their term papers. 

Unfortunately, however, they were not able to 


3.5 Summary o f Student Enrollment Statistics 

An Indication of interest in machine computation may be obtained by con¬ 
sidering the enrollment in each of the subjects which have boon offered during the 
academic year, 1956-1957. 

No. of Students 


Subject 

Title 

Fall 

Spring 

Total 

2.215 

Methods of Engineering Analysis 

17 

0 

17 

6.25 

Machine Aided Analysis 

132 

0 

132 

6.535 

Digital Computer Coding 

0 

112 

112 

6 .538 

Electronic Computational Laboratory 

0 

52 

52 

6.567 

Switching Circuits 

52 

0 

52 

6.568 

Switching Circuits 

0 

48 

48 

15.542 

Management Information Systems 

30 

0 

30 

M371 

Operations Research 

33 

0 

33 

M39 

Methods of Applied Mathematics 

152 

187 

339 

M41L 

Numerical Analysis 

37 

0 

37 

M412 

Numerical Analysis 

0 

18 

18 


TOTALS 

453 

417 

870 
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It is apparent from the above list that the total attendance in these 11 
subjects was 870 students. This doeB not meun however, that there were 87u Indivi¬ 
duals involved, since a number of students took many of the above subjects. Thus 
the actual number of individuals exposed to machine computation is somewhat less. 


K. M. Verzuh 
Assistant Director 
MIT Computation Center 


e 
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1. SYSTEMS ENGINEERING 
WW1 RELIABILITY (18 March - 26 May 1957) 

ptrlod: ^ fOU ° Wln * “ the WWI Computer Reliability for the above-mentioned 

Total Computer Operating Time In Hours 1687 3 

Total Time Lost in Hourt 

40.9 

Percentage Operating Time Usable . 

“7 ,6 

Average Uninterrupted Operating Time Between Failure 
Incidents in Hours 

28.0 

Total Number of Failure Incidents 73 

Failur- incidents per 24 Hour Day ! Q 

Ave- ige Lost Time per Incident In Minutes 31 - 

Average Preventive Maintenance Time Per Day in Hours 2 3 

ANELEX LINE PRINTER 

2 -r r r “ 

printers witHhe^o^tlve^f'cLn^^l' Pt ° Pl0 inve8 “«“‘ e »lgb »P«d line 
quantity of data. The Whirlwind ! people suggested 1 thl^ \ ^° Ut thlS Urge 
was purchased from Anelex Corporation by Division II of ^ “ 

120 characters 56 characters, 

considerations, the printer Is operated at only^SS""•roefT 8 ' 000 ‘° ec * uipmenl 

if all 120 columns per lino are used If per rainute Whirlwind I 

t^speed can increase up to a maximum of ,00 

several days to obtain, "Is^ow 'ready that P r «>vlously took 

notic tape for delayed printout. This results in a ™ v ?!* r ® qUlre<1 to record 0 " »*«- 
as enabling users to have the data available several da^sooner PUt ° r tlne “* * eU 
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2. VISITORS 



a „ _ ‘ -.-saw,., u ,, u an iinuiin.il UlSCUSSlon Ol UU. 

L' T Jt * r '°“ po " ent *- Du ’-ing the past quarter the following * groups totalling 
130 people visited the computer installation. 


April 30 

Cambridge School of Weston 

B 

May 

1 

Evorett Vocational 

High School 

17 

May 

8 

Northeastern University Course 03 21o, Computer 
and Control Components 

25 

May 

13 

MIT Course 16.21 . 

Aircraft Structures 

9 

May 

14 

Brown University 


13 

May 

16 

MIT Course 6.25, 

Machine-Aided Analysis 

35 

May 

22 

Wellesley College 


11 

May 

24 

Representatives of AMEC/Numerical Control 
Subcommittee 01 the Aircraft Industries 
Association 

12 


The procedure of holding Open House at the Laboratory on the first Tuesday 
of each month has continued. A total of 33 people attended the 3 Open House tours 
during the quarter, representing members and friends of the MIT community, Western 
Electric Co., Convair, Northrop Aircraft, Hyde Park High School and Boston Trade School 

Hiring the past quarter, there were also 42 individuals who made brief tours 
of the computer installation at different times. Represented by these Individuals 
were MIT, University of Pittsburgh, Elliott Brothers. American Gas and Electric Co., 
United Aircraft Corp., Harvard University, Michael Baker, Jr., Inc., and Remington 
Hand Corporation. 


3. PUBLICATIONS 


Project Whirlwind technical reports and memoranda are routinely distributed 
to only a restricted group known to have a particular Interest In the Project, and to 
ASTIA (Armed Services Technical Information Agency) Document Service Center, Knott 
Building, Ikivton, Ohio. Requests for copies of individual reports should be made to 
ASTIA. 

The following is a list of memoranda published by the Scientific and 
Engineering Computations Group during the past quarter. 

DCL-161 Generalization of Conjugate Gradient Methods D. N. Arden 1/22/57 
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DCL-165 

SHARE Assembler Operators Notes for Initial 
Use at the MIT Computation Center 

s. F. 

Best 

2/18/57 

DCL-171 

The Direction Cards Used by the Proposed 
Operator Program 

F. C. 

Helwig 

4/18/57 

DCL-I71 -1 

A Proposed Operator Program for the 704 
(Revision of DCL-171) 

F. C. 

Helwig 

5/1/57 

DCL-175 

Summary of PKDS2. an Octal Memory Printout 
Program 

S. F. 

Best 

6/4/57 




PERSONNEL ON THE PROJECT 


COMMITTEE ON MACHINE METHODS OF COMPUTATION AND NUMERICAL ANALYSIS 
Faculty Supervisors: 


Philip M. Morse, Chairman 
SaiF'iel H. Caldwell 
Stephen H. Crandall 
David Durand 
Herman Feshbach 
Jay W. Forrester 
Francis B. Hildebrand 
Chia-Chiao Lin 
J. Francis Kelntjes 

Ex-Officio 


Doan K. Arden 
Fernando J. Corbato 
Frank M. Verxuh 


Physics 

Electrical Engineering 
Mechanical Engineering 
Industrial Management 
Physics 

Industrial Management 

Mathematics 

Mathematics 

Electrical Erglneerlng 


Electrical Engineering 
Computation Center 
Computation Center 


PROJECT WHIRLWIND 

Frank M. Verzuh, Head 
Dean N. Arden 
Sheldon F. Best 
Marion Callaghan 
Frank C. Helwig 
Jack Roseaan 
Arnold Siegel 
Murray Watkins 
Monroe R. Weinstein 
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